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Abstract 

The superfluid properties of the attractive Hubbard model in a Zeeman 
magnetic field, and in the weak coupling regime have been investigated. 
The temperature and magnetic field dependencies of the order parameter 
have been analyzed. Furthermore, the temperature vs. magnetic field 
and temperature vs. spin polarization phase diagrams for the 2D and 3D 
lattices have been obtained. For some parameters a reentrant transition 
has been found. 

1 Introduction 

Spin-polarized superfluidity (in the context of the cold atomic Fermi gases) and 
unconventional superconductivity with a nontrivial Cooper pairing has recently 
been investigated both theoretically and experimentally [1] . 

Due to the presence of a magnetic field, the numbers of particles with spin 
down and spin up are different. This makes the formation of Cooper pairs across 
the spin-split Fermi surface with non-zero total momentum (k f, — k + q |) 
(Fulde and Ferrell [2], and Larkin and Ovchinnikov [3] (FFLO) state) possible. 
The pairing in spin-polarized state is very interesting not only in the context 
of superconductivity, but also in that of trapped unbalanced ultracold Fermi 
atomic gases and color superconductivity in high energy physics [4]- [9]. 

In this paper we analyze the influence of a pure Zeeman effect on the su- 
perfluid characteristics within a lattice fermion (the spin-polarized attractive 
Hubbard) model. For sufficiently high magnetic fields we can observe a change 
in the character of finite temperature transition between the superconducting 
(SC) and the normal state (NS), from the second to the first order. Such a 
phenomenon has been discovered for the first time by Sarma in 1963 [10]. This 
paper consists of three parts. The first part gives analysis of the spin-polarized 
attractive Hubbard model in the Hartree-Fock approximation. The second part 
presents numerical results. In sec. Ill we summarize the discussion. 
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2 Model 

The model Hamiltonian is the attractive Hubbard model (U < 0) in a magnetic 
field [11]: 

ija i i 

where: — c\^c^, = c^Cij, tfj - hopping integral, U - the on-site inter- 
action, fx - the chemical potential. The Zeeman field h can be created by an 
external magnetic field (in (<7/zb/2) units) or by a population imbalance in the 
context of the cold atomic Fermi gases. 

Transforming the Hamiltonian (1) to the reciprocal space, one obtains: 

- ^E(4 T ^-i%)' ( 2 ) 

fc 

where the electron dispersion (with hopping only between the nearest neighbors) 

is eZ = E^iSe* - * = - 2taQ k' 9 fe = T,i=i...d cos ( k i a i) ( here d ^ 2,3 for two- 
and three-dimensional lattice, respectively), ai is the lattice constant in the Z-th 
direction (we set ai = 1 in further considerations). Now, using the Hartree- 
Fock approximation (the pairing only with q = 0), we can obtain the following 
equation for the superconducting order parameter (A = — Efc( c -fcj. c fc-f)) : 



A = -TtV^^ tanh — a + tanh— £± , 3 

iV4"2w £ 2^ 2 2 J w 

where: 

%-H i + ^ T )% + ^ £ + ^ + ^ (4) 

^ T = (-* T ™-h + Ui , (5) 

Ui = y/((-tl-tl)Q £ -ii)* + \A\*, M = M-^> (6) 

= l/k B T, M = n^—ni - spin magnetization (polarization), n a = jj Efc( c j: c k a )^ 
n = n-t + ri| - electron concentration. 



Equation (3) takes into account the spin polarization in the presence of a 
magnetic field and the spin-dependent hopping (t T ^ t^) [12]. 
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The particle number equation takes the form: 

n = 1 - — V — : i-* £ tanh — ^ + tanh — =± . 7 

2N ^ ujj: y 2 2 J w 

The equation for the magnetization is: 




By calculating the partition function in the usual way, we can determine the 
grand canonical potential: 

^ = \vn(2-n)-»+\vMi-¥f (9) 
and also the free energy: F/N = Sl/N + [in. 

3 Results 

We have performed an analysis of the influence of magnetic field on superfluidity 
based on the equations (3 - 9), paying special attention to the behaviour of 
the order parameter and the spin-up and spin-down electron density. In the 
following we set = = t and use t as the unit. 

Fig. 1 shows the temperature vs. magnetic field phase (T — h) diagrams for 
the 2D (a) and the 3D (b) lattices for a fixed chemical potential \i. Both in 2D 
and in 3D, a finite temperature second order phase transition takes place to the 
normal state at sufficiently low values of the magnetic field. With increasing 
magnetic field, the character of the transition between the superconducting and 
the normal state changes from the second order to the first order (the thick 
dotted line), which starts from the tricritical point (TCP). The line of the first 
order transition has been determined numerically from the condition: Q s = Q n 
(where fl n and f2 s denote the grand canonical potential of the normal and the 
superconducting state, respectively). The curve below the first order transition 
line on the phase diagrams (the thin dotted line) is merely the extension of 
the line of the second order transition below TCP. The critical values of the 
magnetic field for a simple cubic lattice are lower than for the square lattice. 
The Hartree term raises the values of the critical magnetic field at T = for the 
first order transition, which exceed the Clogston limit [13] (h c = A /\/2, where 
A is the gap at T = 0, h = 0). 

It is important that the above diagrams were obtained for a fixed chemical 
potential. If the number of particles is fixed instead and n ^ 1, one will obtain 
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Figure 1: Temperature vs. magnetic field phase diagrams for U = —2. (a) 2D 
lattice, ix = -1.358405, (b) 3D lattice, fj, = -f .627404. The arrows indicate the 
tricritical point; the thick dotted line denotes the first order phase transition to 
the normal state. The chemical potentials have been chosen to yield n « 0.75 
at T = and h = 0. 

two critical magnetic fields on the T — h phase diagram [14], [15]. The two 
critical fields define the phase separation region between the superfluid phase 
with the particle density n s and the normal state with the density of particles 
n n . We note that, for a simple cubic lattice and for a weak attraction (Fig. 
lb, U = —2), the differences between the values of these two critical fields are 
negligible and in this case practically only one critical field exists. 

The change in the character of the transition is clearly visible in the behavior 
of the order parameter as a function of temperature, which is depicted in Fig. 
2b. We have the first order BCS to NS transition at T = (Fig. 2a). For low 
magnetic fields, A vanishes continuously with increasing temperature. There 
arise two non-zero solutions for the order parameter, for h = 0.23. The lower 
branch is unstable. We have also found a very interesting behavior of the order 
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Figure 2: Dependence of the order parameter on magnetic field at T = (a) 
and temperature (b) (for three fixed values of the magnetic field), for 2D and 
U = -2, fj, = -1.358405. In Fig 2b, for h = 0.23 the lower branch is unstable. 
For h = 0.24 the vertical dashed line denotes the first order phase transition 
to the normal state and there are two second order transitions in the reentrant 
case (at T = 0.06067 and T = 0.09891). 



parameter for fixed h = 0.24. In this case A vanishes discontinuously at T = 
0.0225 (the solutions with A ^ become energetically unfavorable, i.e. the 
upper branch is metastable and the lower branch is unstable for T > T c .) - thus 
we have the first order transition. However, for T > 0.06067 the superconducting 
solution becomes energetically favorable again (second order transition to the 
SC state) and A vanishes continuously (second order phase transition to the 
NS at T = 0.0998), which is clearly visible in Fig. 3b. Such behavior points 
out that for sufficiently high fields a reentrant transition takes place. Hence, the 
increase in temperature can induce superconductivity. Similar results have been 
obtained for the 3D case (and U = —3). The behavior of the grand canonical 
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Figure 3: Grand canonical potential vs. temperature for h = 0.24 (2D case), 
(a) the first order phase transition to the normal state, (b) the details of the 
reentrant transition: second order transition from NS to SC and from SC to NS 
(inset). The vertical dashed lines mark the phase transition temperatures. 



potential (Fig. 3a) indicates the point of the first order transition from the 
superconducting (Sarma Phase (SP)) to the normal state. 

Figure 4 shows the temperature vs. spin polarization (T—P) phase diagrams 
for the 2D case (a) and the 3D case (b). (P = (nj — nj,)/(nj +nj,)). One can 
distinguish three states for the system under consideration: SP, Phase Separa- 
tion region (PS) and NS. PS terminates at TCP. At T > and P = we have 
the BCS phase. At T ^ 0, A ^ and P ^ 0, the system is in the Sarma phase 
(i.e. superconductivity in the presence of the spin polarization). The transition 
from SP to NS is of second order. Furthermore, we have a transition to the PS 
region (at T = the coexistence of the superconducting (BCS state, A ^ 0, 
P = 0) and the normal state (A = 0, P ^ 0), and for T > the coexistence of 
SP (A ^ 0, P ^ 0) and paramagnetic NS (A = 0, P ^ 0)). Because A vanishes 
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Figure 4: Temperature vs polarization phase diagrams for U = —2. (a) 2D 
lattice, n = -1.358405, (b) 3D lattice, fj, = -1.627404. The dotted lines arc the 
first order phase transition lines, A denotes the gap at T = and P = 0. In 
(b) n w 0.75. 



discontinuously for higher magnetic fields (larger P), the transitions both from 
SP to the PS region and from the PS region to the NS are of first order (the 
dotted lines). Moreover, in the PS region not only the polarizations but also 
the particle densities in the SC and NS are different. 

4 Conclusions 

The influence of the magnetic field on superfluid properties of the attractive 
Hubbard model in a weak coupling regime has been considered. The T — h and 
T — P phase diagrams have been obtained for the two- and three- dimensional 
cases. We can distinguish following phases on the diagrams: the Sarma Phase, 
the Phase Separation region (between the spin-polarized superconducting state 
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(with n. s ) and the normal state (with n n )) and the Normal State. In the presence 
of the magnetic field the densities of states are different for the particles with 
spin down and spin up. Moreover, the magnetic field destroys the superfluidity 
through the paramagnetic effect (or by a population imbalance). One can also 
have a reentrant transition for sufficiently high magnetic fields on the T — h 
and T — P phase diagrams. For 2D system at h = 0, the transition from the 
superconducting to the normal state is of the Kosterlitz-Thouless (KT) type. In 
the case under analysis, the phase fluctuations have been neglected. However, 
for weak attraction the Hartree-Fock critical temperatures may be comparable 
with the KT ones. 

In this paper, we have restricted the analysis to the case of s-wave pairing 
only with q = 0, leaving out the case of the FFLO phase. An important limita- 
tion, as far as the charged superfluid is concerned, is the neglect of the orbital 
effect of a magnetic field and the magnetic field fluctuations [16]. Investigation 
of the competition between the superconducting phases and CDW diagonal or- 
dering (in particular at and close to the half- filled band n = 1) would also be 
an interesting extension of this paper. 
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